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Abstract. We study the statistics of the number of connected components and the 
volume of a random real algebraic hypersurface in RP™ defined by a homogeneous 
polynomial / of degree d in the real Fubmi-Study ensemble. We prove that for the 
expectation of the number of connected components: 

(1) E6 (XBpn(/)) = e(d n ), 
the asymptotic being in d for n fixed. 

We do not restrict ourselves to the random homogeneous case and we consider more 
generally random polynomials belonging to a window of eigenspaces of the Laplacian 
on the sphere S n , proving that the same asymptotic holds. As for the volume properties 
we prove that: 

(2) EVol(Z RP n(/)) = 6(d). 

Both equations and |2| exhibit expectation of maximal order in light of Milnor's 
bound bo(Zup»(ff) < 0(d") and the bound Vol(Z KP «. (/)) < 0(d). 



1. Introduction 

Hilbert's sixteenth problem, in its general form, asks for the study of the maximal 
number and the possible arrangements of the components of a nonsingular real algebraic 
hypersurface of degree d in WP n . It turns out that even the case of curves in MP 2 is an 
extremely subtle problem and essentially nothing is known for n > 3. The analogous 
question for a hypersurface in CP n is trivial, as for a generic choice of the polynomial 
defining it, its topology is determined (this is the statement that the cohomology of a 
complete intersection in the complex projective space is determined by its multidegree 
and complete intersections are generic). 

Thus, if from one hand no technique exists to study the arrangement of the compo- 
nents of a real algebraic hypersurface, on the other hand one still would like to formulate 
"typical" statements: a way for doing that is to replace the word generic with random. 

Let us start by considering the case of hypersurfaces in MP 1 , i.e. zeros of ordinary 
univariate polynomials. In fact here we consider the homogenization of a univariate 
polynomial together with its projective roots rather than the classical polynomial itself, 
but it is clear that the two objects produce the same statistic. Over the complex numbers 
statements like: 

"a generic complex polynomial of degree d has d roots" 

make perfect sense. On the other hand the only property we can say to be generic for 
the real roots of a real polynomial is the one of being distinct. 

l 
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Proposing a random viewpoint raises the obvious question: what is meant by random? 
We first look to the seminal work [H] of Kac, who proved that a random univariate 
polynomial of degree d whose coefficients are i.i.n. has asymptotically 2/7rlog(i real 
roots. Even if it might seem natural, the distribution Kac considered is not invariant by 
the action of projective transformations over the extended real line. To produce such an 
invariant distribution we proceed as follows. We first extend our space of polynomials 
by homogenization; this leads us to consider projective zeros rather than just on the 
real line. Thus, we consider the vector space W\ d of real homogeneous polynomials of 
degree d in two variables. We endow this vector space with a Euclidean structure by 
considering the L 2 -norm of the restriction of these polynomials to the circle S . In an 
equivalent way we can define for / and g in W\d their scalar product to be: 

(f,9)= I f(x)g(x)e-W\ 2 dx, 



where x is the variable in M 2 . Notice that since / and g are homogeneous of the same 
degree, the previous integral equals the integral of fg restricted to the circle S 1 (up 
to a constant factor). Using this scalar product we can define a centered Gaussian 
probability distribution on by letting for a measurable subset A: 

probability of A = / e ~ df, 

v l,d J A 

where v\& = (2-7r)( d+1 )/ 2 is the constant that makes the previous a probability distribu- 
tion and / is treated as the integration variable that runs over W\d- 

For us a random polynomial will simply be a polynomial in with the above 

distribution. 

Having thus replaced generic with random, we can make a positive statement in the 
real setting: 

"a real random polynomial of degree d has ^ ^ ^ roots". 

We provide now another equivalent definition of a random polynomial. We fix a ba- 
sis {/i, . . . , fd+i} for Wi t d which is orthonormal for the above scalar product and we 
consider a random function of the form: 

d+1 

/ = with £i~#(o,i). 

1=1 

What we get in this way is again a random polynomial in the above sense. One way to 
produce such a basis is by considering trigonometric polynomials, i.e. functions on the 
circle satisfying: 

A s if = -l 2 f, l en. 

We denote by H\ i the vector space of all the solutions of the previous equation: these 
must be linear combinations of sm{W) and cos(/#). Equivalently we can extend such 
functions by homogeneity to all M 2 and we get harmonic polynomials of degree I. What 
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is interesting for us is that in this fashion we can write the vector space W\ t d as: 

Wi, d = \\x\\ d - l H hl , 

(2-Z62N 

where again x is the variable in M 2 and the direct sum turns out to be orthogonal with 
respect to the above scalar product. 

In particular the above orthonormal basis can be obtained by collecting together the 
orthonormal bases coming from each H\i, d — I £ 2N. 

If we repeat the same construction for the scalar product 

</,<7>C = / f(z)g(z)e-^ 2 dz, 
Jc 2 

then the corresponding real random polynomials are called Kostlan distributed (or 
Bombieri-Weyl) . The vector space is the same, but in this second case the Euclidean 
structure comes by restricting to the real polynomials the Hermitian structure given 
on the complex ones. The expectation of the number of real roots of a real Kostlan 
polynomial of degree d is Vd. Already from this comparison we see that the structure 
of the zeros of a random polynomial is richer than the one of a Kostlan, which in turn 
is richer than for the naive ensemble. 

Remark 1. The two numbers d(d+2)/3 and d are called the parameters of these distribu- 
tions. More generally each Gaussian distribution on W\ t d that is invariant by projective 
transformation has its parameter 5 which is defined to be the square of the expected 
number of zeros it produces. 

Our interest is in what happens when we increase the number of variables. We let 
W n ^ be the vector space of homogeneous polynomials of degree d in n + 1 variables, 
and consider their zero sets as hypersurfaces in WP n . If in the univariate case as a 
"measure" of our set we could use the number of its points, here we have at least two 
possible approaches to generalize this idea. 

The first one is by considering the number of points as a zero-dimensional volume. In 
this case an upper bound for the "size" of our set is provided b}[j] 

(3) Vo1(2hp»(/)) < G(d). 

Remark 2. The previous formula comes from an application of the Integral Geometry 
formula and Wirtinger's inequality (see |2Uj). In fact if / is regular the previous volume 
is bounded by d ■ Vol(MP n ^ 1 )/Vol(CP n - 1 ) and since d = Vol(Z CP ™(/)) this also gives a 
nice way to compare the real zero locus to the complex one (as we did in the univariate 
case using the Fundamental Theorem of Algebra). 

The second approach is more subtle and regards the number of points as a topological 
property (rather than a metric one), namely as the number of connected components. It 
is interesting to notice that there is a quantitative bound also for this number, usually 
referred to as Milnor's bound (see [19]): 

(4) b Q {Z W n{f)) < Q(d n ). 

^Here the notation f(d) = Q(g(d)) means that there exist two constants ci,C2 such that for d large 
enough cig(d) < f(d) < C2g(rf). 
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Remark 3. The original paper [19] of Milnor does not contain the bound in this form, 
since a d n+1 term appears. We will give a proof of Q in the Appendix for the smooth 
case (which is the case we are interested in, since the set of nonregular polynomials has 
measure zero in W n j), though the same bound holds true even for singular hypersurfaces. 

As we already pointed out, for a generic / in W n d both the volume and the Betti 
numbers are constant for its complex zero locus. This is certainly not the case for the 
real zero locus, as for example it might even be empty. To study the expectation of 
these quantities over the reals we make W n ^ into a Gaussian space. We first define the 
scalar product of two elements / and g in it to be: 



Using this scalar product W Ut d is turned into a probability space by defining for each of 
its measurable subsets A: 



where now u n4 = (2vr) dim ( w/ «^)/2 > 

Following [23], we call W n d with the resulting distribution of probability the real 
Fubini-Study ensemble. 

The expectation of the volume of the zero locus of a random polynomial of degree d 
from the real Fubini-Study ensemble was computed by P. Biirgisser in [TJ: 



In fact in the paper [7J this result follows from the more striking computation of the 
expected curvature polynomial of Zkp«(/) for any projectively invariant distribution. 
Such a result is formulated in terms of the above mentioned parameter 5, which in this 
case is the square of the average number of real zeros of a restriction to an MP 1 . 

In particular, if we define the Kostlan distribution on W n <j by means of the scalar 
product^] J* c „ + i f(z)g(z)e~^ z ^ 2 dz, the corresponding parameter is d and the result we get 
for the expectation of the volume is 



(the general result for a distribution with parameter 5 is <5 1/ ' 2 Vol(IRP n ~ 1 )). 

In the case n = 2 and restricting the distribution of probability to #2,d ; i-e. consid- 
ering only harmonic polynomials rather than homogeneous ones, the statistics of the 
number of connected components was studied by Nazarov and Sodin. In the break- 
through paper [22J the authors prove that there exists a constant c > such that for a 
random spherical harmonic / of degree d: 




probability of A 






EVol(Z KP n(/)) = VdVoKMP™- 1 ) 




Moreover, they showed that bo concentrates exponentially about its mean. 



'The corresponding norm is usually called the Bombieri norm. 
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In the visionary handwritten letter [23] P. Sarnak claims that using the same ideas as 
in [22] it is possible to extend this result to the case of real homogeneous polynomials in 
the real Fubini-Study ensemble. The interesting part of these statements (besides the 
concentration) is the lower bound estimate: an upper bound of the same order is given 
by Milnor's one (or in the spherical harmonic case by Courant's nodal domain theorem). 
The statistics of the number of connected components of a Kostlan distributed curve 
was studied by Gayet and Welschinger: in the paper [TJ] the authors prove that such a 
real curve of degree d has on average 0(d) components. The same authors prove in |12] 
that the expectation of each Betti number of a Kostlan hypersurface in MP n is bounded 
by 0(d n / 2 ), generalizing in some sense the univariate case. 

On a different perspective the first author in [16] . using techniques from Random 
Matrix Theory, has studied intersection of one or two Kostlan distributed quadrics, the 
asymptotic being this time in the number n of variables; for such random algebraic sets 
the expectation of the sum of the Betti numbers behaves asymptotically as n. 

As we did before, it is possible to give an equivalent definition of the real Fubini-Study 
structure on W n d using an L 2 (S' n )-orthonormal basis. In fact we can again decompose 
our vector space into the orthogonal sum: 

Wn,d= M\ d ~ l H n ,U 
d-le2N 

where now H n i is the vector space of spherical harmonics of degree I, i.e. those functions 
on the sphere S n satisfying: 

A S nf = -l(l + n-l)f. 

Inspired by the ideas of [22] and [23J we study a more general class of random polyno- 
mials that contains as special cases both the homogeneous and the harmonic ones. To 
define such a class, for every a G [0, 1] and d > let us consider the following index set: 

[ad, d] = {natural numbers I such that [ad] <l<d and d — I is even}. 

With this definition the a- window of eigenspaces we are going to consider is the following 
orthogonal direct sum: 

H [a d,d\ = H n<u a G [0,1]. 

l£[ad,d] 

We endow H[ ad d ] with the probability distribution induced by the real Fubini-Study 
ensemble. 

Remark 4. There is no unified convention for the names of these ensembles. In fact the 
name real Fubini-Study should be distunguished from complex Fubini-Study which is 
another name for what we called Kostlan distributed. 

In this paper we will only deal with the real Fubini-Study ensemble and its restriction 
to -ff [c^d] • For this reason, since no confusion can arise, a random polynomial for us will 
simply mean a polynomial from this ensemble. 

In particular for a = we get random homogeneous polynomials of degree d and for 
a = 1 we get random spherical harmonics of degree d: 

H [o,d] = w n,d and H[d,d] = H n,d- 
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Remark 5. In fact we do not need to restrict to the case the a-window contains only 
the eigenspaces H n j with d — I even. Our choice is motivated by the fact that we want 
to avoid redundancy in order to get random polynomials for a = 0; one can consider as 
well the full window. 

We denote the dimension of H n i by d(n, I) and the one of Hy ad ^ by D a ; notice that in 

the case a/lwe have D a = Q(d n ). If we consider an orthonormal basis {Xi}t=i ^ or 
H nt i, we can decompose a random / in i?r a <i,<i] into a sum of random spherical harmonics 
of different degrees: 

d(n,l) 

f= E with zt-mD- 1 ). 

l£[ad,d] i=l 

The introduction of the scaling coefficient D~ l is absolutely not necessary, but it has 
the advantage that E||/||^ 2 (gn) = 1 which simplifies the notation in the proofs. 

As for the expectation of metric properties of the zero locus of a random / in H\ ad ( j\, 
we prove the following theorem. 

Theorem 1. Let a G [0, 1] and f be a random polynomial in Ht a( £ )( n. Then: 

EVol(Zs»(/)) = 0(d). 

In particular taking expectation on both sides of inequality ^ turns it into an equality 
(the volume on the projective space and the one on the sphere being related by a constant 
that depends only on n). 

We notice that the statement for a = 1 was proved by Berard in [5] and for a = by 
Biirgisser [7| as stated above in ([5]). In fact it is possible to write the exact value of this 
expectation using an analog of what we called the parameter of the distribution, defined 
in this case as 5 = ^jy^ Ylie[add] 1(1 + n — l)d(n,l); the expectation of the volume is then 
^VolOS" 1 - 1 ). The same statement holds for the projective zero locus, as it is double 
covered by the spherical one. 

The statistics of the number of connected components is more subtle and requires 
more work, drawing on the ideas of [22]. The statement we will prove is the following. 

Theorem 2. Let a G [0, 1], and f be a random polynomial in Ht a( i t( n. Then: 

Eb (z S n(f)) = e(d n ). 

In analogy with the case of the volume, this theorem says that if we take expectation 
on both sides of Milnor's bound Q we turn it into an equality (the number of connected 
components on the projective space and those on the sphere being related by a factor 
of two). 

Notice that the case n = 2 and a = 1 gives the spherical harmonic case proved by 
Nazarov and Sodin in [22J. The case n = 2 and a = gives Sarnak's claim [23J on curves 
from the real Fubini-Study ensemble. 

Remark 6. After this paper was written, the authors have been informed by F. Nazarov 
and M. Sodin that they have an alternative approach to this problem, that works in 
more general settings; a summary of their ideas is available at [24j . 
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The authors have also been informed by P. Sarnak that together with I. Wigman they 
are able to provide asymptotics for the expected number of ovals of a random plane curve 
for the naive model (the defining coefficients for the monomials are independent normals 
with mean zero and variance one), the Kostlan model and the real Fubini-Study one. 

1.1. A random version of Hilbert's sixteenth problem. It is interesting to discuss 
the case a = and n = 2, i.e. random real algebraic curves. The first part of Hilbert's 
sixteenth problem asks to study the configuration of the components of a real algebraic 
curve of degree d in the plane MP 2 . The nullhomotopic components of such curves 
are called ovals; each one of them separates MP 2 in two components: one of them is 
homeomorphic to a disk (and is called the interior of the oval) and the other to a 
Moebius band. In the case d is even all components are ovals, but if d is odd there is 
a component that is not an oval (it is homotopic to a projective line and generates the 
fundamental group of MP 2 ). Harnack's theorem states that such curves cannot have 
more than ( rf ~ 1 K rf ~ 2 ) _|_ \ components, but the possible arrangements of their ovals is a 
far more complicated problem; we refer the reader to the beautiful survey |26| on the 
subject. 

In fact the problem itself of constructing curves with many ovals is a difficult one. 
Gayet and Welschinger's theorem [lOj on the exponential rarefaction of maximal curves 
confirms such difficulty: if we write the equation of a curve putting some coefficients in 
front of the standard monomials, it is exponentially difficult to get one with many ovals. 

This "approach" of building curves depends on the choice of the basis we want to 
use for if instead of the monomial basis (which is orthogonal for the Bombieri 

scalar product, inducing the Kostlan distribution) we use the spherical harmonics the 
result improves dramatically. In fact now, because of Theorem [2] for a = and n = 2 
the expectation of the number of components of the curve we obtain is a fraction of 
the maximal one. The price we pay is that we are no longer allowed to use the simple 
monomial basis and we need an explicit expression for the polynomials representing 
spherical harmonics, such an expression being quite complicated. 

A completely analogous discussion holds for real hypersurfaces: now the fact that 
what we get applying Theorem [2] is on average a fraction of the maximal number of 
components is due to Milnor's bound. 

The method employed here (based on the ideas of Nazarov and Sodin) suggests that 
more generally one could ask for the expected arrangement of such ovals. In fact essen- 
tially the proof relies on the following observation: given any point x on the sphere and 
a small disk around it of radius 0(d _1 ) there is a positive probability (independent of d) 
of finding an oval inside such small disk (and such oval loops around x). Unfortunately 
this probability, even if constant, is still very little. A numerical computation of the 
expectation of the number of components for n = 2 was performed by M. Nastasescu in 
|21| : as Sarnak writes in his letter "a random curve is 4% Harnack". Thus, the random 
version of Hilbert's Sixteenth Problem, to study the expected arrangement of a random 
curve, is still open and despite its difficulty it seems to be very fascinating. 

One of the main difficulties in this framework is to define a statistic (beyond the 
number of components) that captures some features of the arrangement. We suggest 
one that we call the energy of the curve C. 
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First we say that one oval C\ is contained into another one C2 if it lies in the interior 
of Ci. Using containment we can give a partial ordering on ovals and construct a forest 
(a collection of trees) representing their arrangement. The vertices of the trees represent 
the ovals, with the empty ovals represented by leaves, and an oval not contained in any 
other is represented by a root. We define an energy that is additive for disjoint unions 
of trees, and multiplicative when appending trees. So the energy of a forest is the sum 
of the energies of its trees. We define the energy of a tree with only one isolated vertex 
as 2, and of a larger tree inductively as twice the sum of the energies of all subtrees that 
result from removing its root. This description determines the energy of any forest (for 
example, see Figure [I]). 




Figure 1. A hypothetical arrangement of ovals for C and the corre- 
sponding forest. The energy is h(C) = 2 • 2 • (2 • 2 + 2) + 2 • (2 + 2 • 2) = 36. 

By Bezout's Theorem the length of each tree cannot be more than |_f J and if such 
a tree exists then no other trees can appear; the corresponding curve is said to be 
maximally nested and its energy is h(C) = 2 d l 2 . On a somehow opposite direction one 
can consider Harnack curves, having h(C) ~ 5c? 2 (for d large enough). Such curves 
have the maximal number of components and satisfy some special condition on the 
arrangement of their ovals; we refer the reader to [18J for the exact definition and to [15] 
for surprising connections with planar dimers. 

Thus, using this notion, one open problem that we can state on the expected ar- 
rangement is to describe the asymptotic behavior as d goes to infinity of the energy of a 
random curve. The pictures we get of a random curve suggest that the energy has the 
same order as a Harnack curve, so we conjecture that: 

lim Y**m = 2 

d->oo log d 

(the lower bound follows from Theorem [2]). 

1.2. Plan of the paper. In Section 2, we specialize to the case n = 1. This helps to 
illustrate the main techniques used in the paper and provides a proof for the value of 
the parameter 5 used Section 6.2. Theorem [T] is proved in Section3 and Theorem [2] in 
Section 4. Section 5 is devoted to the proof of the main technical lemmas for Theorem [2] 
In Section 6 we discuss some examples (more on random curves and random surfaces). 
In the Appendix we provide the proof of Milnor's bound, we review a construction of 
hyperspherical harmonics, as well as some of their properties. 
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2. A TOY MODEL PROOF: THE UNIVARIATE CASE 

To illustrate the main ideas, we first consider a random polynomial in the special case 
n = 1. As mentioned in the introduction, in this case both the zero dimensional volume 
and the number of connected components coincide. For simplicity of notation we will 
denote their expectation by: 

Ed = Ejnumber of zeroes of a real random univariate polynomial}. 

In view of the generalization to higher dimensions, we will give two proofs that E& = 
Q(d). The first one is based on techniques coming from integral geometry and will 
produce an exact result; this proof will generalize to any number of variables producing 
again an exact result. The second is based on the barrier method introduced by Nazarov 
and Sodin; it will produce only an approximate result but has the advantage that it 
generalizes to higher dimensions. 

We notice that in the case n = 1, a real random polynomial is just a trigonometric 
polynomial with either only even or only odd terms (depending on whether d is even or 
odd). 

2.1. Integral geometry proof. The first proof is obtained following the ideas of [9]. 
Recall that if we have functions {fi, ■ ■ ■ , fk} defined on the real line we can consider the 
random function: 

k 

/ = with &~#(o,i). 

1=1 

Then if the coefficients are independently distributed, the expected number of zeros of 
/ on the interval [a, b] is given by the formula: 

1 f b 

(6) Ejnumber of zeroes of / on [a, b}} = — ||7(t)||di, 

7T J a 

where 7 is the projection of the curve c : 1 1->- . . . , fk{t)) on the unit sphere S k ~ 1 

(we refer the reader to the beautiful paper [9j for this and related topics). We note that 
^ is proved using a special case of the kinematic formula from integral geometry. In 
our case [a, b] = [0, 2tt] and if d is even, 

% 1 sr~^ ( sm(2kt) , cos(2/rf) \ . , „. 
f(t) = £ ^= + V a k \_ ; + b k ^-L , £0, a fc A ~ N(0, 1 . 

The curve c is given by: 

c(t) = —j= ^— j=, sin(2t), cos(2t), . . . , sin(id), cos(td) 
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and its norm in M fc is constant and equals \ /^jr- In particular we have: 



The norm of 7 is 

||7|| = ||(0,2cos(2*),-2sin(2*),...,dcos(td) ) -dsiii(td))[| 
which easily simplifies to: 



d+V 



3 

A similar calculation in the case d is odd gives the exact same result. We are now in 
the position to use formula Q and we obtain that the expected number of zeros of / on 

[0, 2tt] is 2W ^ rf ( rf + 2 l \ Thus, the number Ed of expected zeros on the real projective line 
is: 



/«(<* + 2) 



3 

Notice that this is an exact result and gives the asymptotic Ed = Q(d). 

2.2. Barrier method proof. In this section we only prove that Ed > cd for d large 
enough; the inequality Ed < d follows from the Fundamental Theorem of Algebra. The 
proof is centered around the function described in the following Claim whose proof is 
deferred until the end. 

Claim (Existence of the barrier). There exist numerical constants p and c\ such that, 
for each sufficiently large d and each x G S 1 there is a homogeneous polynomial b x of 
degree d and I? -norm one satisfying: 

b x (x) > Cl Vd and b x {e ip,{d+l) ■ x) < —c\ Vd. 

Once a point x in S 1 is fixed, we can decompose / = ^ob x + f x where £0 is a Gaussian 
random variable with variance ^rf, and f x is in the orthogonal complement to b x in 
W\ d- We choose a new Gaussian random variable £0 independent of £0 and of f x with 
variance ^jj-j-, and set: 

f± = ±6A + fx- 

By construction the random functions /+ and /_ have the same distribution as /. We 
notice that we can write: 

f = Hab x + l(U + f-). 
Next we apply the following Claim to /+ and /_. 

Claim. There exists a constant C2 such that K\f(x)\ < C2 for every x in S 1 . 

We also defer the proof of this Claim in favor of first seeing how it is used. Let us recall 
that Markov's inequality for a positive random variable \X\ reads IP{|^| > c} < c _1 E|X|. 
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If we apply such inequality to \X\ = because of the previous Claim we can find 

a constant Co such that: 

P{|/±(x)|>cn}<l/5 and P{|/ ± (e-/( d+1 )x)|>C }<l/5, 

where p is provided by the Existence of the barrier Claim. 

The event that |/±(x)| and \f±(e lp ^ d+1 ^x)\ are each simultaneously at most Co is an 
intersection of four events each of which is complementary to the one of the above. Thus 
the probability of their intersection satisfies: 

P (n) = 1 - P (U) > 1 - 4 • 1/5 = 1/5. 

We note that these four events are not mutually independent, but this is not needed in 
the line above. Now consider the event Q(x) that f(x) > cq and f(e tp ^ d+l ^x) < —cq. 
If Q(x) occurs then / has a zerd^] on S 1 between x and e lp ^ d+1 ^x. The event Cl(x) 
happens provided that: (i) ^ov^ > 2d[ cq and (ii) f± evaluated at at x and e tp ^ d+1 ^x 
are each bounded in absolute value by cq. Note that (i) and (ii) are independent, and 
we have just checked above that (ii) happens with probability at least 1/5. Moreover, 
PjCo 1 /^ > 2cJ~ 1 co} > C3, for some constant C3, since ^oyfd has variance ^q^r = 0(1). 

It remains to choose Q(d) disjoint intervals in S 1 each of length p/d. Each of them 
contains a zero of / with probability at least C3/5 and hence: 

E d >c-d, 

as desired. 

First Claim. We will choose our barrier in such a way that the barrier b x centered at 
some point x G S l is obtained just by precomposing fr(o,i) with a rotation that moves 
the point (0, 1) to x, so we only need to define &(n,i) which we do by normalizing the 
function: 

tt 1 i0\ tt 1 a ■ a\ sin((<i + 1)6>) 
U d {e w ) = U d (cos6,sm6) = — . - . 

smt) 

Notice that U d is a polynomial of degree d, since the denominator divides the numerator. 
(Ud is a projectivized Chebychev polynomial of the second kind). 
We have: 

U d (0, 1) = d, 

say by de l'Hopital's rule with 9 — > 0, and for all d sufficiently large, we also have 

u d (V 5 ^) < --d. 

Thus it suffices to show that the norm of U d is on the order of \fd as d — > 00, for then 
we will obtain the desired properties for 6( 10 ) (the normalization of U d ). 
It is easy to show that ||C/d|| 2 = 2nd; thus the function 

h^) = ^f d 



In this one dimensional case the existence of a zero is essentially the intermediate-value property. 
In higher dimensions this event will be generalized by considering the boundary of a disk centered at x; 
in the current situation we only need to use one of the two boundary points of an interval centered at x. 
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1 „„J „ _ 3n 



has the desired properties with c\ = y= and p 



2V^k a ^ H 2(3+1) • 

□ 

Second Claim. Let us choose a new orthonormal basis go, g%, .., g<i so that all basis ele- 
ments except go vanish at x. To see that it is possible to do this, first notice that the 
subspace V 1 - of homogeneous polynomials of degree d vanishing at x has codimension 
one inside Wi t d- Thus, to get the desired orthonormal basis, we may first orthonormal- 
ize a basis for V ± , then take its orthogonal complement V which consists of just one 
function. Once normalized, this function gives the basis element go. 
Writing our random polynomial in terms of this basis, we have 

1/0)1 = \Co9o(x) + ^2^k9k(x)\ = |£o| • \9o(x)\- 

k 

Using again the standard basis fk, the fact that go is a normalized trigonometric poly- 
nomial of degree d means that 



bo 0)1 



^2 A kfk(x) 



< Yl \ A k\ < E W^d+l = Vd+l, 



where we have also used the triangle inequality and the discrete Cauchy-Schwarz in- 
equality. Combining this with the above, 



|/(x)| < Iv^+Wol, 

which implies that E|/(x)| is bounded by a constant (the Gaussian random variable 
\/d + 1 • £o has variance of constant order, and therefore so has the expectation of its 
modulus) . 

□ 

3. Proof of Theorem Q] 

We generalize the proof given in [5] for the harmonic case to our a-window of 
eigenspaces. Here the main idea is to use a basis of the eigenfunctions of the Lapla- 
cian to define a kind of Veronese embedding, transforming the geometry of zeros of such 
functions on S n to the geometry of the intersection of the image of S n with hyperplanes. 
What is nice about using harmonic functions is that such a Veronese embedding turns 
out to be a dilation, i.e. dilates the Riemannian metric by a constant factor. In the 
language of integral geometry such embeddings are called moment maps (generalizing 
the idea of the moment curve 7 that we used in the univariate case in Section 2.1). 

Let us denote by d(n, I) the dimension of H n j as above (thus in the case n = 2 we have 

d(l,2) = 21 + 1 and in the general case d{n,l) = e(/ n " 1 )). For each (n,l) let {Yi)f={ l) 
be an orthonormal basis for H n j and consider the map: 

: S n R d(n ' l) 

whose components are (Y"i, . . . , Ym u n). Such a map has the following properties: (i) 
it is an immersion and its image is a submanifold of the sphere in M d ( n ' 1 ^ of radius 
R = (d(n, /)/Vol(5' n )) 1 / 2 ; (ii) it is a dilation of ratio XR 2 /n where A = 1(1 + n - 1) it is 
the eigenvalue for the eigenspace H n i of A^n . In other words the length of the image of 
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a unit vector in S n under the differential of A n> j is \J\/nR. Property (i) follows from 
the addition theorem for spherical harmonics: 

d(n,l) 

Z l (x,y)= Yi(x)Yi(y), x,y€S n . 
i=i 

Here Z\ denotes the (reproducing kernel) zonal harmonic of degree I, for which Zi(x, x) = 
Vol n (s n ) • Thus, evaluating the addition formula at y = x we get: 

d(n, I) 



Yol(S r 



|A„,/(x)|| 2 



which shows the image is contained in the sphere. Property (ii) is proved in [6] and goes 
under the name of nice imbeddings property. 

For a given a G [0, 1) we consider thus the map: 

A : S n -> R Da = R d{n ' l) 

l£[ad,d] 

whose components are the A n> / normalized by the factor {D a /Vo\{S n )) 1 / 2 , i.e. 



fVol(S r 



OL 



1/2 



A=(— ) (A ad ,...,Ai,...,A d ), d-leven. 



With this choice of the normalization the addition theorem for spherical harmonics 
implies that the image of A is contained into the unit sphere in R Da . We claim now that 
A is a dilation of ratio 5, as defined in the Introduction. To prove this let us consider a 
unit vector v £ TS n and let us compute the square of the norm of its image under the 
differential of A: 

le[ad,d] a le[ad,d] a 

where in the second equality we have used the fact that A n j is a dilation of ratio 
X(n,l)R 2 /n. We apply now the integral geometry formula (see [9]): 

(7) A *n _ n ( Vol({a}^nA(5"))da = Vol(A(5")) VolLS! ' ^ 



VoliS *- 1 ) J S D a -i VL J v " K v " Vol(5") ' 

Property (ii) above implies that the map A is a covering space (in fact either a double 
cover, or a diffeomorphism) ; let us say it is /c-sheeted. We claim now that the left hand 
side of Q equals exactly k5 l ~ n - 1 ^ 2 ¥No\{Z sn {f)). In fact if a = (a x , . . . , a Da ), then 
{a} 1 - n A(S n ) equals the image under A of the nodal line of the function Yld=i a iYii 
where {1^}^ is the basis for H\ a d,d] obtained by collecting together the elements from 
the various orthonormal basis of H n j. Since A is a covering, the volume of such image 
(which is n — 1 dimensional) equals kd^ 1 ^ 1 ^ 2 times the original volume; moreover since 
our Gaussian distribution is uniform on the unit sphere in Hy ad d ^ , then the above integral 
computes expectation and the conclusion follows. 
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As for the right hand side of ([7]) again the fact that A is a /c-sheeted covering and 
a dilation of ratio 5, implies Vol(A(S' n )) = kd n / 2 \o\(S n ). Plugging all this into (m) we 
finally get: 

EVol(Z s »(/)) = 5 1/2 \o\{S n ~ l ). 

4. Proof of Theorem [2] 

Our proof is deeply influenced by the ideas and the exposition of [22j. In fact for 
a = 1 our claim is the direct generalization of the lower bound estimate in [22J; this case 
is a simple extension of the proof in [22] and we devote a remark to it at the end of the 
proof. 

Notice that the upper bound for our expectation is provided by Milnor's bound Q, 
discussed in the Appendix, and we only have to prove the lower bound. 

Thus we fix^Jfor the rest of the proof an a G [0, 1). To start with, for x a point on the 
sphere S n we consider the following event: 

(8) fi(x,r) = {/(z) > and /| aD ( x , r ) < 0}. 

If 0(z, r) occurs, then {/ = 0} has a component inside D(x, r) (and this component loops 
around x). If we prove that for r = @(d~ 1 ) there exists a constant a± > independent 
of x and d such that 

F{tt(x,r)} > oi, 

then we can cover the sphere S n with Q(d n ) disjoint such disks and get the statement 
(each one of these disks contributing at least a\ to the expectation in the statement). 
Specifically, the radius r = r(d, n) we are going to use from now on is given by: 

_ %n 
T ~ 2d + n- 1' 

where 1/2 is the first point of minimum of the Bessel function Jo and y n for n > 2 is the 

first point to the right of zero where the Bessel function J n-2 reaches a minimum (in 

2 

this case J n-2 already has a minimum at zero, see 1251 ) . 

2 

The following Lemma provides the existence of the so called barrier function; we defer 
its proof until the next section. 

Lemma 3 (Existence of a barrier). For every a G [0, 1) there exists a constant c\ = 
c\{a) > such that for every d > and every x in S n there exists a homogeneous 
polynomial B x of degree d and I? -norm one satisfying: 

B x {x)> Cl d n / 2 and B x \ dDM <- Cl d n / 2 . 

Using B x we can decompose H^ ad d ] = span{i? x } © sp&nlBx} 1 - , thus getting the 
following decomposition for /: 

f = ^oB x + f ± with Co- N^D- 1 ). 



In particular every constant in this proof is allowed to depend on a, but not on d. 
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We let now £o be a random variable distributed as £o but independent of it and we 
define f± = ±£oB x + /-*-. Notice that both / and f± have the same distribution. The 
introduction of these new random polynomials allows us to write: 

(9) f = CoB x + ^(f + + f-) 

and to split our problem into the study of the behavior of B x and f± separately. In 
fact the event £l(x, r) happens provided that for some constant 02 > the two following 
events both happen: 

1) E(x,r) = {CoB x (x) > 2a 2 and £oB x \ dD (x,r) < -2a 2 }; 

2) G[x,r) = {||/±(a?)[| < a 2 and \\f±\dD(x,r) II < 02}- 

To check that E{x,r) n G(x,r) implies Q(x,r) we simply substitute the inequalities 
defining E(x,r) and G(x,r) in ^ evaluating respectively at x and at dD(x,r). 
Now by definition the two events E(x, r) and G(x, r) are independent and thus: 

F{E(x, r) n G(x, r)} = ¥{E(x, r)}F{G(x, r)}. 

Because of this it suffices to show that the probability of each one of them is bounded 
from below by a positive constant that does not depend on d and this will provide a 
bound from below for the probability of Q(x, r) which is independent of d. 

The fact that the probability of E(x, r) is bounded from below independently on 
d immediately follows from Lemma 3| since D a = Q{d n ) and B x {x) = ®{d n / 2 ) the 
random variable £,qB x (x) is Gaussian with mean zero and variance B(l); similarly for 
the boundary. Thus for 03 small enough the probability of it being bigger than 03 is 
uniformly bounded from below. 

It remains to study G(x,r); to this extent notice that this event is given by the 
intersection of the four events G\ )2 = {||/±(^)|| < 0-2} and (^3,4 = {||/±|aD(a;,r) II — a 2}- 
For these we do not need independence as: 



f{ n G t } = 1 - p{ ( n G t y} = 1 - p{ g > 1 - net] 

i i i i 
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and it is enough to prove that the probability of the complement of each one of them is 
small. 

To bound the probability of G\ and G\ we argue as follows. First notice that since 
f± are distributed as / it is enough to prove the corresponding bound for a random 
/. We consider in H\ ad ( j\ the hyperplane V 1 - of those functions vanishing at x together 
with its orthogonal complement V. Such orthogonal complement must be spanned by 
a norm-one function F of the form: 

„2 



(10) 



F 



< 1, 



Fo+ E mYl with E y i 

l£[ad,d] le[ad,d] 

where Yi £ H n j is the normalized zonal harmonic centered at x (notice in particular 
that this normalization produced the r.h.s. inequality in the above line). In fact in 
each space H n i the orthogonal complement of the zonal Yi centered at x vanishes at 
this point, and we immediately get (10). Such zonal function is obtained by composing 
the function in equation (19) with a rotation that moves the north-pole of S n to x; 
moreover by setting 9 = in equation (19) we see that it has order 0(^ n_1 )' 2 ) at x. 
Using the above orthogonal decomposition H\ a( i t d] = V ® V 1 - , we can write the random 
/ as ^qF + F 1 - with F 1 - a random function in V 1 - . Thus we get: 



(11) 



E|/(aO| <E 



£o E \ymx)\ 

l£[ad,d] 



< E 



^0 



l£[ad,d] 



1/2 



< 03- 



We have used ( 10 ) and the discrete Cauchy-Schwarz inequality in the second inequality; 
for the last inequality we have used the fact that Yiix) 2 = ©(Z n_1 ), which implies 

1/2 



s{d) 



l£[ad,d] 



6(o m / 2 ) 



Since the variance of £ has order 0(d n ), when we scale it by s(d) we make its variance 
of constant order, which gives (11). The statement on the probability of {|/(^)| > ^2} 
being smaller immediately follows from Markov's inequality. 

We explain now the argument for the remaining two cases. Again since f± is dis- 
tributed as /, the probability of each G| 4 is exactly the probability of the following 
event: 

{[|/Uz>(a!,r)ll > «2 for a random / G H [ad4] }. 
Hence we consider the positive random variable: 

I XI = max I f I. 

dD(x,r) 

Bounding the expectation of \X\ is more delicate and we devote to it a separate lemma, 
whose proof is postponed to the next section. 

Lemma 4. For every a 6 [0, 1] there exists a constant C2 = 02(a) such that 

oU(x,r) 

Applying Markov's inequality to \X\ combined with the previous Lemma gives a 
bound for the probability of G| 4 , concluding the proof. 
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Remark 7. This proof does not apply when a = 1, but in that case one can use a rather 
direct generalization to n + 1 variables of the proof given in |22j . Namely, their Claim 3.1 
providing a bound on the expected max over a small disk is still true in more variables 
and can be proved in the same way, and the barrier function can again be chosen to be 
the normalized zonal harmonic composed with a rotation. 



5. Proofs of Lemmas 

5.1. Proof of Lemma [3], We work out the details for n = 2 first. We will construct 
B x for x the north pole and define it for the other points by considering the composition 
with rotations (such rotations preserve all the properties of the statement). 

Set y = y2 (the first minimum point of the Bessel function Jo) and let YJ° be the 
normalized zonal harmonic of degree I centered at the North Nole. Using the Legendre 
polynomial Pi of degree I we can write: 

Y?{0,<l>) = P,(cos0)V2Z + l, 

where (9, 0) are the standard polar coordinates on S 2 . We recall from |25j the following 
asymptotic for Pi (Hilb's asymptotic): 

where Jo is the Bessel function of the first kind and order zero. The error Ri{6) is always 
less than Q 1 l 2 0{l~z). The function Jo satisfies the following properties: Jo(0) = 1, it is 
decreasing in [0,y] where y is a local minimum with Jo(y) < 0. We call z £ (0,y) the 
first zero of Jo. Let be the maximum of the two values a and z/y; we define B x as: 

^2lG[nd,d] 

B x — 



a/ Cardf^ci, d] 

Since each Y® has norm one and since different zonals are pairwise orthogonal, we 
immediately see that the L 2 -norm of B x is one, as required from the statement. Let us 
now estimate the value of B x at x. Using the integral estimate for the numerator, which 
provides a bound of the form 6(d 3 / 2 ), and the fact that ^Cardf/wi, d] = @(d 1 ^ 2 ), for 
some constant a± > independent of d we have: 

B x {x) = > aid. 

a/ Card[//d, d] 

Recall that the number r was defined such that 2rf 2 l " 1 r = y, i.e. we set: 

2y 

r = — : . 

2d+l 

We estimate now the value of B x on the boundary of a disk of radius r. Notice that it 
suffices to show the bound for the numerator in the definition of B x has order 0(d 3 / 2 ), 
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since we already know the denominator has order Q(d 1 ^ 2 ). Because of Hilb's asymptotic 
we have: 

B{d 1 l 2 )B x {r) = Y, \[^ Jo (^V ir ) VW + l+ E Ri(r)V2l + l. 

le[fid,d\ * V / l£[fj,d,d\ 

v v ' v v ' 

L(r) R(r) 

We estimate the error R(r) first. Notice that the summands Ri(r)\/2l + 1 are all 
bounded by r 1 / 2 0(/ _1 ) so that: 



R(r) < r l l 2 \J2 r j = rl/2 0{logd) = O 



log d 
Vd 



Next we bound L(r) as follows. First notice that since fi = max{a, z/y}, then for 
every I in the sum Jq{^^t) < 0. We fix now a 7 G (fx, 1) and notice that since each 
term is negative and Jo is decreasing in (7?/, y), then: 

/ \ V 2 

L(r)< E — Miy)V2i + i. 

^-^ \ sm r I 

le[yd,d] v 7 

Moreover using the expansion (r/ sinr) 1 / 2 = 1 + Q(d~ 1 ) we get: 

L(r) < J ( 7 y)(i + e(^ 1 ))0(d 3 / 2 ) < - C2 d 3 / 2 , 

where we have used the integral estimate for the sum of the square roots. Recalling that 

L(r) = B x (r)Q(d 1 ^ 2 ), this concludes the proof for n = 2. 

For the case n > 2 we argue as follows. First recall that we set y n to be the first point 

after zero where the Bessel function J n-2 reaches a minimum. The function J n-2 has the 

2 2 
following properties: it vanishes at zero and we denote by x n its first point of maximum 

(it is positive at this point); it is decreasing on the interval [x n ,y n ] and negative at 

y n . We let z n be its only zero in [x n ,y n ] and fi n be max{a, z n /y n }. If Y\ S H n< i is the 

normalized zonal harmonic centered at x, we define the barrier as: 

Yl,l£\ii n d,d} Yi 



a/ C&rd[fj, n d, d] 



Since the zonals are pairwise orthogonal, the L 2 -norm of B x is one. As already discussed 
in the proof of Theorem|2l the order at x of each zonal Y\ is B(^ n_1 ^ 2 ) as shown in ( 18 ), 



and thus the order of their sum at x is 9(d( n+1 )/ 2 ). Moreover since ^/Cavd[fi n d, d] = 
®(d l l 2 )i then we immediately get that for some positive constant c±: 

B x {x) > Cl d n ' 2 . 

The proof of the bound for the value of B x on dD(x, r) involves the same steps as in the 
case n = 2, this time using a generalization of Hilb's asymptotic involving the Bessel 

function Jn-2 . We need an explicit expression for Yj first: 

2 

( n ~2 n — 2 \ 

Yi{9) = c(n,l)P; 2 ' 2 j (cos 9), 
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where P l 2 ' 2 is a Jacobi polynomial and c(n,l) has order 0(l 1 / 2 ), see (20) in the 
Appendix. 

Using the asymptotic provided in ([25] , Theorem 8.21.12) we have: 



(n^n^) / I) l)\ r / \ 1 / 2 



P, v 2 ' 2 y (cos0) = ( sin--cos-j \ h ( n > l ){—^) Jn=2(N9) + R l {9)\, 



where N = ^f=^, /i(n,Z) = 9(1) is given by (plk, and the error term Ri(9) is always 
less than 9 1 / 2 0{l^). 

Since by definition r = we have (sin | • cos §)^~ = 0(rf~2~), and plugging all 

this into the definition of B x , we have that B x {r) equals: 



n — 2 , 



C / \ 1/2 

(-)(,/-,,(-)(,/--) £ (_!_] J ZL _,(r(2/ + »- l)/2) + /?,(,-■; 



smr 



For the error term we get a bound of the form 0{d 2 ). As for the principal term, we 

fix as above a 7 n G (/z n , 1) and notice that since each summand is negative and J n-2 is 

2 

decreasing in (j n yn,y n ), then: 

^(r)<e(^) V Q(l ll2 )J^i( ln yn)<-cid n l 2 
» — • 2 

le[f n d,d] 

where in the last inequality we have used the integral estimate again. 

5.2. Proof of Lemma |4[ Since the distribution is invariant by rotation, it is enough 
to prove the statement in the case the point x equals the north-pole (0,0, ..,0, 1). We 
need an explicit expression for a real orthonormal basis of H^ add y To this end, we use 
the inductive construction of an orthonormal basis for H n j, the space of hyperspherical 
harmonics in n+1 variables of degree I, given in terms of orthonormal bases {Yj((f>)}jei m 
for each space H n _i^ m with m = 0,1,.., I, where I m is an index set of size \I m \ = 
dim H n -i )Tn = d(n — 1, m). Namely, the basis elements as reviewed in the Appendix are: 

YT(e,<f>) = Nr(smerp(sf +m \ C o S e)Y0). 

In the previous formula pj^ is the Gegenbauer polynomial [251 Formula (4.7.1)], 9 is the 
angle from the north-pole, the multi-angle G gives the remaining coordinates, 

and the normalization constant N™ is stated explicitly in (16) (see also [2j Chapter 3] 
for an exposition of this inductive construction). 

Using this basis we can write the random / in H[ add j as: 

/= E EEtfiHM) 

111 

= E E E W-T +m W)>m 

le[ad,d] m=0j£l m 
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with each $ distributed as iV(0, D^ 1 ), with D~ l = <d(d~ n ). In particular, the restriction 
of / to dD(x,r) is given by: 

/M) = E E { E W-T +m) (co S r)}(sinrry,(0), 

where the index set L m := {/ G [ad, d] 1 1 > m} arises from having reversed the order of 
summation. 

Notice that the term in the large brackets is a Gaussian random variable itself: it 
is a linear combination of independent random variables, each one of mean zero and 
variance D^ 1 . We denote such a random variable by ^; in other words we define: 

(12) i j = £ ^p/_X +m) (cosO, 3 e Im- 

By the addition formula for independent Gaussians we get that $ is distributed as a 
Gaussian with mean zero and variance a(m,d) 2 where: 

(13) a(m,d) 2 < £ {ATp/_X +m W)} 2 - 

l=m 

Remark 8. Most importantly, for each m the Gaussians £ J are identically distributed 
across all j within the index set L m . This follows from the fact that in the sum (12) the 
coefficient of depends on m but not j. 

Using this new notation we can rewrite: 

d d 
f( r ^) = E E e{^r) m Y 3 ^) = ^(sinr) m £ gYj (</>). 

m=0j£l m m=0 jelm 

Using the triangle inequality, we bound the expectation of \X\ = maxgj}( Xj1 .) \ f\ as: 



d 

E\X\ < V(sinr) m E max 

771=0 



E ^ 



By the previous remark the function: 

f(<p) = J2 i jY M 

is a random spherical harmonic on S'™ -1 with independent identically distributed Gauss- 
ian coefficients. This allows us to apply the following Claim, whose proof is provided in 
the Appendix. 

Claim (A basic estimate for random spherical harmonics:). 

E max \F(4>)\ < a(m,d) • d(n — l,m), 
(pes™- 1 

where recall that d(n — l,m) is the dimension of the space of spherical harmonics of 
degree m in n variables. 
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Returning to E|X|, we have: 



(14) 



E|X| < (sin r) m a(m,d) -d(n-l,m). 



m=0 



It remains to estimate a(m,d). We note that the Gegenbauer polynomials satisfy 
Formula (7.33.1)]: 



max 

xe[-i,i] 



P 



l—m 



(X) 



Using this to estimate (13) we get: 



mm, 



l=m 



I + m + n — 2 
I — m 



I + m + n — 2 
I — m 



Inspecting the constant (N™) 2 (see formula (16)), we notice that we can bound part of 
it by a constant: 

n+2m+l ■j 



' n+2m 



)(n + 2m-l) 



(in fact this number goes to zero as m goes to infinity); in particular we bound each 
term in the above sum by a constant times 

2 



(n + 2m - 2)1(2/ + n - - m)\ (I + m + n - 2 

I — m 



which simplifies to 
This finally gives: 



(Z + m + n-2)! 
(2l + n-l) 



I + m + n — 2 
I — m 



a(m,d) 2 < ^-^(2Z + n-l) 



D, 



l=m 



I + m + n — 2 
I — m 



< 



2ai 



< 



< a 



D a ■ (2m + n-2)l < 
2ai 

D a ■ (2m + n-2)! 
(2d) 2m 



^2(1 + m + n - 2) 2m+n - 1 
(d + m + n- 2) 2m+n 



2 (ml) 2 



where a<i is a constant that does not depend on d. 

In particular a(m, d) < a 2 ■ Recalling that r = 2 d+n-i — d ( we have set y n = y 
for simplicity of notation) and using (sinr) m < r m < (|)"\ we have in (14): 

•' (2y) r " 



m=0 



< (I) ■o(m,d)-d{n-l,m)<a 2 Y j 



m=0 



m\ 



■ d(n — 1, m). 
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We estimate d(n— 1, m) < a^m n 2 < ci4 m , so that using the definition of the exponential 
function we obtain: 

E |X| < a 5 y < ase 2 ** < c 2 . 

^— ' ml 

m=0 

This concludes the proof. 

6. Examples 

6.1. More on the arrangements of random curves. Besides the number of ovals 
of a curve in the projective plane, there are other interesting invariants associated to it 
that are worth studying (e.g. the energy we defined in the Introduction). 

One classical invariant is the number of empty ovals; let us denote it by uq. If the 
curve C has even degree d = 2k, then Arnold's inequalities (see [1] section 6) imply: 

MC) > b (C)-(k-l)(k-2). 

In particular if a curve has the maximal number of ovals, then vq{C) > k 2 . The numerical 
results obtained by M. Nastasescu in |21j have shown that the constant in the bound 
for Kbo is very small: 

Ebo(C) < ^d 2 , 
25 

and combining this with Arnold's bound would only give ¥.vq(C) > 0. On the other 
hand, as stated in the Introduction, the barrier method can be used to give information 
also on the expected arrangement of the components of the random hypersurface. In fact 
the event 0(a;, r) as defined above Q implies the existence of an empty oval inside the 
disk D(x,r). Thus, placing on the sphere 0(d 2 ) disjoint such disks gives the following: 

Average number of empty ovals: Ei/o(C) = 0(<i 2 ). 

Remark 9. We notice that every / 6 W n 2k defines a function on the real projective 
space WP n and for a generic / the gradient of this function does not vanish on its zero 
locus. In particular if C is a component of ZjRp™(/), then MP n \C has two components; 
if q : S n — > WP n is the quotient map, we define the interior of C to be the component 
of MP n \C whose preimage under q is not connected. As we did for curves, we say that 
a component C\ of Z^p n {f) contains another C2 if C2 is contained in the interior of G\\ 
an empty component is one whose interior is connected. 

Using this notation the above corollary can be stated for hypersurfaces as well and 
provides the existence of 0(d n ) empty components. 

Similarly we can define the energy of an hypersurface as we did for curves: additive 
on disjoint unions and multiplicative on nestings. Instead of always taking the value 
2, the "seed" we take for the energy of an empty component is the sum of its Betti 
numbers (in the case of an empty oval, this equals 2). 

As for higher moments of bo, Nazarov ad Sodin have proved that for a = 1 (random 
spherical harmonics) it exponentially concentrates about its mean and Sarnak has in- 
formed the authors that the same result holds for a = (homogeneous polynomials). 
Thus the typical curve has 0(<i 2 ) components. The authors expect the same results to 
hold for any n as well. 
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On a different direction Gayet and Welschinger have studied random curves in the 
Kostlan ensemble (the complex Fubini-Study one). In [10] they proved that maximal 
curves become exponentially rare in their degree. Specifically for any sequence {a(d) > 
l}deN* ° ne can consider P(a(d)), the probability that a Kostlan curve has more than 
(d-i)(d-2) _|_ ^ _ components, and there are positive constants C and D such that: 

P{a(d)) < Cd 6 e~^. 

For this random model they have also proved that the expectation of the number of com- 
ponents is less than 0(d), conjecturing that this was the right asymptotic (as confirmed 
by Sarnak). 

Finally, for curves there is another random model introduced by the first author in 
|16j that has not been discussed so far and has relation with Random Matrix Theory. 
Recall that once we fix three symmetric matrices Aq, A±,A 2 of order d, we can consider 
the following homogeneous polynomial: 

f(xo,x 1 ,x 2 ) = det(x A + xiAi + x 2 A 2 ). 

If we let the symmetric matrices to be independent random matrices in the Gaussian 
Orthogonal Ensemble, f becomes a random polynomial of degree d; its distribution is 
very different from the Fubini-Study ones. Despite this, every real algebraic curve arises 
as the zero locus of one such polynomial and we can define random curves in this way; we 
call them random determinantal curves. For example, using this description, maximally 
nested curves correspond to matrices for which there exists a positive definite linear 
combination cqAq + c\A\ + c 2 A 2 and since the probability of the positive definite cone 
becomes exponentially small in d, they are very unlikely (notice that using the above 
definition they also have very high energy). 

Remark 10. Each symmetric matrix A of order d also defines a quadratic form a : M. d — > 
M. by the following identity: 

a(x) = (Ax, x) for all x £ R d . 

It is interesting that under this correspondence A is in the Gaussian Orthogonal En- 
semble if and only if a is Kostlan distributed. Using this property the first author has 
computed the expectation of the sum of the Betti numbers of the intersection of one or 
two Kostlan distributed quadrics in IRP^" 1 , showing that it behaves asymptotically as 
d (see [16]). 

Moreover there is a kind of duality between the curve above defined and the common 
zero locus of ao,a\,a 2 , as discussed in [T7] . This duality allows to interchange the 
topology of the random intersection of three quadrics and the one of corresponding 
random determinantal curve, i.e. to study one using the other. The authors plan to 
discuss this ensemble and related properties in a forthcoming paper. 

6.2. Betti numbers of random surfaces. Consider a random polynomial / £ W^d 
and its zero locus Y = Z s -<,(f) on the sphere S 3 . Such Y is a random surface and 
using the above results, combined with the computation of the expectation of its Euler 
characteristic, we can provide an asymptotic for each of its Betti numbers. Applying 
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Theorem 1 from J7J with s = 1 and n = 3 we get: 

E X (% P 3(/)) = -^(^-3) ! 

where 5 is the parameter of the distribution, which we computed to be d ( d ^ , Since 
Y double covers Z KP 3(/), its Euler characteristic is twice the one of Z KP ?,{f). Thus, 
plugging the value of the parameter in the above formula we get: 

E X (Y)~-^= = e(d 3 ). 

We notice now that since Y is the zero locus of a globally defined function on S 3 , then 
it is orientable (a nowhere vanishing normal is given by grad(/)|y). We can thus apply 
Poincare duality and get the equality bo(Y) = b 2 (Y). Using the definition of the Euler 
characteristic we get: 

Eb (Y) = Eb 2 (Y) = G(d 3 ), 

and 

E6i(Y) = E(6 Q (K) + b 2 (Y)) - E X (Y) = G(d 3 ). 

We notice that in the case / has even degree then it defines a function on MP 3 as well, 
thus Z KP 3(/) is orientable and the same asymptotic holds; in the case d is odd, the zero 
locus on the sphere is still orientable, but the one on the projective space might have 
nonorientable components; in particular Poincare duality holds only for 7L 2 coefficients. 
In this case we get similar asymptotic using Betti numbers with 7L 2 coefficients. 

7. Appendix 

7.1. Milnor's bound. As we said in the Introduction, we only need to prove the bound 
for a regular / (the corresponding zero locus is smooth); the set of nonregular / is a 
proper algebraic subset of W n d, hence with measure (and probability) zero. In fact 
this proof can be modified to obtain a bound of the same order also for nonregular 
polynomials, using techniques from semialgebraic geometry, but this is far beyond our 
scope. 

To start with, notice that it suffices to give a bound for the sum of the Betti num- 
bers, since this will bound in particular 6o(-^RP n (/))j moreover the Universal coefficients 
Theorem implies that the sum of the Betti numbers with integer coefficient is bounded 
by the sum of the Betti numbers with Z2 coefficients (see [13f ) : 

&o(2bp»(/)) < b{Z UP n{f),Z) < b(Z RP n(f),Z 2 ). 

Now we use the so called Smith's inequality, as presented in [26]. This is the state- 
ment that for a topological space X with an involution r : X — >• X, the inequality 
6(Fix(r),Z2) < b(X,Z 2 ) holds. In particular, since / is a real polynomial, the complex 
conjugation is an involution on ^cp™(/) an d Smith's inequality reads: 

6(Z KP n(/),Z 2 ) <&(2cp»(/),Z 2 ). 
We use now the fact that Zcp u (/) is a nonsingular hypersurface in CP™ of degree 
d. This property has two consequences: first, its cohomology has no torsion and 
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&(Zcp«(/),Z2) = 6(Zcp™(/),Z) (Lemma 3.1 of [8]); second, such cohomology is com- 
pletely determined by d and n, and in particular: 

6(Zcpn(/),Z 2 ) = 6(Zcp»(/) > Z) 

= (d ~ 1)n+1 d ~ ( ~ 1)n+1 + n + = e(cP). 

The explicit formula in the above line is a consequence of the adjunction formula, which 
provides the Euler characteristic of Z<cp™(/) (see Exercise (3.7) in [8]); the formula for 
the sum of the Betti numbers immediately follows from Poincare duality. Finally, since 
Zs n (f) double covers Zrp"(/), the same bound holds for: 

b (Z S n(f)) < 2b (Z MP n(f)). 

It is interesting to notice that in the case / is harmonic, this bound can be obtained 
using Courant's Nodal Domain Theorem. In fact this theorem implies that the number 
of nodal domains of / on S n , i.e. bo(S n \Zsn(f)), is bounded by Q(d n ) and the long 
exact sequence for the pair (S n , S n \Zs"-(f)) provides the bound (the reader is referred to 
[1] for these effective techniques from homological algebra). This proof using Courant's 
theorem applies even when the zero set of / is not regular, but it doesn't work for a 
linear combination of harmonic functions. 

7.2. Hyperspherical harmonics. In this section, we review a few properties of the 
hyperspherical harmonics, including any relevant properties of the Gegenbauer polyno- 
mials. 

First recall that the Gegenbauer polynomials are defined using the Jacobi or- 
thogonal polynomials pj> a '^ with a = (3 = A — 1/2 |25l Formula (4.7.1)]. Namely, 

pW(x)= 5 (A,0P/ A - lAA - 1/2) (x), 

where 

T(A + l/2) r(Z + 2A) 



r(2A) r(/ + A + i/2)' 



As mentioned in Section 5.2, an orthonormal basis for H n [ can be built up inductively. 
Namely, start with an orthonormal basis for 

Hn—l,m- 

m£[0,l] 

For each such basis element, we obtain a basis element for H n \. 

This entails that the dimension of the space H n i is the same as the dimension of the 
space © m£ [Q n i? n _i jm . It is instructive to check this using the known formulas for each, 
that is, to verify directly that 

. (n + l\ ( n + I — 2\ (n + m— 1\ fn + m — 3 

„ )-{ n „_! )-{ „_! 

which follows from cancellation in the sum along with the recurrence relation for binomial 
coefficients. 
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We follow [2j Chapter 3] in order to state the details resulting from the inductive 
construction of an orthonormal basis for H n j (restricted to the sphere S n ) of degree I. 
Suppose an orthonormal basis {Yj(cj>)}j & i m for i? n _i ]m has already been constructed for 
each m = 0,1, ..,1, where I m is an index set of size \I m \ = dim// n _i jm = d(n — l,m). 
We can then obtain a spherical harmonic in n + 1 variables (restricted to the sphere S n ) 
of degree / as follows 

(15) Yr(e,J>) = Nr(sm9) m P^ +m \™sO)YM- 

In the previous formula P^ is the Gegenbauer polynomial, 9 is the angle from the north- 
pole, the multi-angle 4> £ S'™ -1 gives the remaining coordinates, and the normalization 
constant NJ 71 is given by: 



(16) N t 



r( n+2 2 m+1 )(n + 2m - 2)!(2Z + n - - m)\ 
\' V^r(^f^)(n + 2m- + m + n - 2)\ 



T m _ . / ■*■ V 2 _ 

iv 

2 

Writing these functions for all j £ I m and for all m = 0, 1, .., d generates an orthonormal 
basis for H n ^. 

In or der to produce the normalized zonal harmonic Yi(9) that was heavily utilized in 

2 



r(==i) 



n—l 



Section |5.l| take m = in the formula (15). So the polynomial Yj((p) = ^ 
is constant and: 

(17) Y t (9) = Y?{9) = ^p} 2 j (cos0), 

where w n _i is the (n — l)-dimensional volume of <S n_1 and 

/ r(^)(n-2)!(2/ + n-l)/! 
i y V5fT(5)(n-l)(i + n-2)l ' 

The Gegenbauer polynomials satisfy [25, Formula (4.7.3)]: 

4*) (1)= («+»-*). 

Together with (17) this implies that 

(is) y,(o) = e(/("- 1 )/ 2 ). 

/ n — 2 re — 2 ^ 

Writing this in terms of the Jacobi polynomial P, 2 ' 2 , 

/ n — 2 re — 2 \ 

(19) Y0) = c(n,l)P; 2 ' 2 '(cos 9), 
where 

(20) c(n,0 = ^< 7 f r ^ 1 ,z)=e(/V2 ) . 

OJn-l V 2 y 

Theorem 8.21.12 in [25J states an asymptotic for the Jacobi polynomials. In partic- 
ular this applies to the polynomial P^ A ^ 2,A 1 ^ 2 \ and taking A = r ^ L , we obtain the 
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asymptotic that was used in Section 5.1 



(n^n^) / I) () \ C / ft \ 1 / 2 



P; 2 ' 2 ; (cos6) = ^sin--cos-J |/i(n, I) ^ — ) J^{N9) + 12,(0) }>, 
where 

iV= 2/ + n-l 



2 

2 — n 

, , ,x {2l + n- 1\~ T(l + n/2) 
(21) fc(n, I) = f J fi = ^)' 

and the error term -R;(#) is always less than O^^OilT^) (see the remark in [25] immedi- 
ately following Theorem 8.21.12). 

7.3. Proof of the basic Claim stated in the proof ofLemma|4) Using the notation 
of Section 5.2, let us recall the Claim: 

E max \F\ < a(m, d) ■ din — 1, m). 

gn-l 

The proof uses a standard reproducing kernel argument. Let Z m (r],4>) be the zonal 
harmonic in n variables with pole at the point <p G S™^ 1 . This is not the normalized 
zonal, it is the reproducing kernel as defined in [3, p. 94]. Namely, Z m (7],(j)) acts as 
a reproducing kernel for the whole space of homogeneous harmonics of degree m. In 
particular, 



F(<f>) = / F(n)Z m (r),<t>)d V . 
The Cauchy-Schwarz inequality yields 

(22) \F{(t>)\ 2 < d(n- l,m) / F(r]) 2 dr], 



where we have used f S n-i Zd(Vi ^) 2 dr] = d(n — l,m) [3, p. 95, Proposition 5.27]. 

Now, taking the maximum over the sphere, then taking expectations over both sides 
of the inequality (22), we have: 

E max |F| 2 < d(n - l,m)E / F(r]) 2 dri 
= d(n-l,m) 

< d(n — l,m) 2 a(m, d) 2 . 

The inequality in the Claim now follows from the general inequality E|y| < y^Ejyp for 
a random variable Y (it is an application of Cauchy-Schwarz again). 
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